COLOR AND SPIN IN QUARKONIUM PRODUCTION" 



ERIC BRAATEN 
Department of Physics, Ohio State University 
Columbus, OH 43210, USA 
E-mail: braaten@pacific.mps.ohio-state.edu 

I describe the NRQCD factorization approach to the inclusive production of heavy 
quarkonium, contrasting it with the color-singlet and color evaporation models. 
These approaches differ dramatically in their assumptions about the roles played 
by color and spin in the production process. They also differ dramatically in their 
predictions for the production of charmonium at large transverse momentum. 



1 Introduction 

The production of heavy quarkonium in high energy collisions is important, 
because there are several quarkonium states with clean experimental signatures 
through their decays into lepton pairs. Measurements of the production rate 
of these states can be used to test our understanding both of heavy-quark 
production in high energy collisions and of the formation of bound states from 
heavy quark pairs. The NRQCD factorization approach, which is based on the 
use of an effective held theory called nonrelativistic QCD to exploit the large 
mass of the heavy quark, has led to dramatic progress in our understanding of 
inclusive quarkonium production. Below, I describe the basic ideas underlying 
the NRQCD factorization approach, emphasizing its implications for the role 
played by color and spin in inclusive production. 



2 Simple Production Models 

Before 1995, most of the effort to understand the inclusive production of char- 
monium was carried out within either the color-singlet model or the color evap- 
oration model. Both of these models have roots that go back to 1975 shortly 
after the discovery of charmonium. These two models make diametrically op- 
posite assumptions about the roles played by color and spin in the production 
process. 
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2.1 Color-singlet Model 



In the color-singlet model, charmonium states are interpreted as nonrelativistic 
bound states of a c and c interacting through a confining potential. The J '/ip, 
for example, is identified as a bound state of a cc pair in a color-singlet 3 Si 
state. I denote such a state by cci( 3 Si), where the subscript specifies the 
color state (1 for singlet, 8 for octet) and the argument specifies the angular- 
momentum quantum numbers. The \cj is identified as a cci( 3 Pj) bound state. 
The only role of gluons in the color-singlet model is to generate the potential 
that binds the c and c. 

The color-singlet model provides a prescription for calculating not only the 
inclusive production rate for a quarkonium state, but also its inclusive decay 
rate into light hadrons and its decay rate into leptons or photons. For example, 
the formula for the decay rate of the J/ ip into light hadrons in the color-singlet 
model is 



where T is proportional to the annihilation rate at threshold of a cc pair in 
the state cci( 3 Si). The last factor in (|l]) is the square of the wavefunction at 
the origin, which gives the probability for the cc pair in the J/ip to be close 
enough to annihilate. 

In the color-singlet model, it is assumed that a cc pair that is produced 
in a high energy collision will bind to form a given charmonium state only if 
the cc pair is created in a color-singlet state with angular-momentum quantum 
numbers that match those of the bound state. For example, to form a J/ip, 
the cc pair must be produced in a color-singlet 3 <Si state. The formula for the 
inclusive production cross section for the J/i/j in the color-singlet model is 



where a is proportional to the production rate for a cc pair at threshold in 
the state cci( 3 S\). The last factor in (0) is the probability that a pointlike cc 
pair in the state cci( 3 5'i) will bind to form a J/tp, which is again given by the 
square of the wavefunction at the origin. 

The color-singlet model is an extremely economical phenomenonological 
framework for calculating quarkonium production. The factors T and a in 
(1) and (2) can be computed using QCD perturbation theory in terms of the 
charm quark mass m c and the running coupling constant a s (m c ). The only 
other phenomenological parameters are a single wavefunction factor for each 
spin multiplet: Rj/^(Q) for the S-wave states J/ip and rj c , R' x (0) for the P- 
wave states \c0, Xci, Xc2, and h c , etc. Moreover, these wavefunction factors 



T{J/i>) = f(cc!( 3 Si)) 
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can be determined from experimental measurements of the decay rates. Thus, 
the color-singlet model predicts the inclusive production cross sections for all 
quarkonium states in terms of a single phenomenological parameter m c . 

The color-singlet model gives definite predictions for the dependence of the 
cross section on the polarization of the charmonium state. It also predicts that 
ratios a(H)/a(J/ifj) of the cross sections of quarkonium states with different 
J quantum numbers should vary dramatically from process to process, be- 
cause of angular momentum selection rules. Unfortunately, these predictions 
have proved to be wrong. Dramatic variations in the ratios <r{H) / a(J /ip) 
have not been seen experimentally and there is little if any evidence for spin 
asymmetries in the cross sections. 

The most dramatic failure of the color-singlet model came in 1995, when 
the CDF collaboration measured the production cross sections for charmonium 
states at the Tevatron pp collider. They used a vertex detector to separate 
prompt production of charmonium from production through the decay of b 
quarks. They found that the cross sections for the direct production of J/ip 
and ip' at large transverse momentum were larger than the predictions of the 
color-singlet model by about a factor of 30u This dramatic discrepancy marked 
the final demise of the color-singlet model. May it rest in peace. 

2.2 Color Evaporation Model 

In the color evaporation model, the probability that a cc pair produced in a 
high energy collision will bind to form a specific charmonum state is assumed 
to be almost completely independent of the color and spin state of the cc 
pair. The basic assumption is that the exchange and emission of soft gluons 
destroys any correlations between the color and spin state of the cc pair when 
it is created and the quantum numbers of the final cc bound state. In the color 
evaporation model, the formula for the inclusive cross section for producing a 
J I ip is conventionally written in the form 

cr(J/V>) = d{cc :4ml<s< Am 2 D ) fj^, (3) 

where a is the cross section for producing a cc pair with invariant mass below 
the DD threshold. This cross section is summed over all the color and spin 
states of the cc pair. The factor /m in (3) is a phenomenological parameter 
that gives the fraction of cc pairs below the DD threshold that form a J/ip. 

The color evaporation model is much less ambitious than the color-singlet 
model. It makes no attempt to relate production cross sections to annihilation 
decays. However it is still a fairly economical phenomenological framework 
for calculating quarkonium production. The cross section a in (3) can be cal- 
culated using perturbative QCD as a function of m c ,m£i, and a s (m c ). The 
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choice of Amr D as the upper endpoint for the integration over s is completely 
arbitrary, so we should regard s max — Am 2 D as a phenomenological parame- 
ter. The remaining phenomenological parameters are a fraction fu for each 
quarkonium state H. Thus the color evaporation model gives predictions for 
inclusive production cross sections in all high energy processes in terms of m c , 
•Smax, and a fraction fjj for each quarkonium state. 

The color evaporation model gives simple predictions for the ratios of cross 
sections for different quarkonium states and for the dependence of the cross 
sections on the polarization of the quarkonium state. It predicts that ratios of 
cross sections are independent of the production process: 



The model also predicts that cross sections are independent of the polarization 
of the charmonium state. For example, the fraction of J/ip's that are trans- 
versely polarized is predicted to be 2/3 for any production process. Thus the 
color evaporation model can be ruled out by experimental measurements of a 
nonzero spin asymmetry for the J/tp in any process or of a variation of the 
Xc-to-J/ip ratio in different processes. 

The color evaporation model is not yet dead. It is basically a phenomeno- 
logical model, so it can only be killed by experimental data. I believe that its 
demise is only a matter of time until sufficiently accurate experimental data is 
available. The model is based on far too simplistic a picture of the effects of 
soft gluons on color and spin in quarkonium production. 



3 The NRQCD Factorization Approach 

The NRQCD factorization approach to quarkonium production was developed 
by Bodwin, Braaten, and Lepage in 19950 As in the color evaporation model, 
a cc pair in any color and spin state has a nonzero probability of binding to 
form a given charmonium state. As in the color-singlet model, that probability 
depends strongly on the color and angular-momentum state. Thus the roles 
of spin and color in this approach are somewhere intermediate between their 
roles in the color-singlet and color evaporation models. However, in contrast 
to these models, the NRQCD factorization approach is based firmly on QCD. 



3.1 Quarkonium as a Nonrelativistic Bound State 

The NRQCD factorization approach to annihilation decays and inclusive pro- 
duction is based on the fact that heavy quarkonium is a nonrelativistic bound 
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state. Its properties are determined by QCD in terms of essentially only 2 
parameters: the QCD coupling constant a s and the heavy quark mass M. 
However, because it is a nonrelativistic bound state, there are a number of dif- 
ferent energy scales that play an important role in quarkonium physics. These 
scales include 

• M, the mass of the heavy quark, which sets the scale for the mass of the 
bound state, 

• Mv, the typical momentum of the heavy quark in quarkonium, which 
sets the scale for the size of a quarkonium state, 

• Mv 2 , the typical kinetic energy of the heavy quark, which is also the 
scale for splittings between radial excitations and between orbital angular 
momentum excitations, 

• Aqcd, the scale of nonperturbative effects associated with light quarks 
and gluons. 

The scales Mv and Mv 2 owe their existence to the fact that a small parameter 
v is generated dynamically in a nonrelativistic bound state by the balance 
between kinetic energy and potential energy. If the interactions between the 
quark and antiquark are described by a potential V(R), the balance can be 
expressed in the form 

Mv 2 ~ V(R) fori? - \/{Mv). (5) 

If the mass M is large enough, the potential is essentially Coulombic: V(R) ~ 
a s (R)/R. The condition (5) then reduces to v ~ a s (Mv). If the mass M is 
not so large, the potential may be dominated by the linear confinement region: 
V{R) ~ k 2 R. In this case, the condition (5) gives v ~ (k/M) 2 / 3 . In cither 
case, we get a dynamically-generated small parameter v. 

In Table [t| I list the energy scales M, Mv, and Mv 2 for charmonium and 
bottomonium and also what they would be for toponium if the top quark did 
not decay too rapidly to form bound states. For the scale M, I have taken 
half the mass of the lowest quarkonium state. For the scale Mv 2 , I have 
taken the average of the energy splittings for the first radial excitation and 
the first orbital angular momentum excitation. The scale listed for Mv is just 
the geometric mean of M and Mv 2 . By dividing Mv 2 by M, we find that 
v 2 is about 1/3 for charmonium, 1/10 for bottomonium, and about 1/100 for 
toponium. While 1/3 is not tiny, it is small enough that we can hope to use it 
as an expansion parameter for relativistic corrections. 
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Table 1: Momentum scales for quarkonium systems 





cc bb 


tt 


M 


1.5 GeV 4.7 GeV 


180 GeV 


Mv 


0.9 GeV 1.5 GeV 


16 GeV 


Mv 2 


0.5 GeV 0.5 GeV 


1.5 GeV 



In Table |], I list the running coupling constant of QCD at each of the 
momentum scales for charmonium, bottomonium, and toponium. For charmo- 
nium and for bottomonium, the scale Mv 2 is definitely in the strong-coupling 
region. At the scale M, a s is small enough that we should be able to use 
perturbation theory to calculate the effects of that scale. 

3.2 Nonrelativistic QCD 

The large mass of the charm and bottom quarks gives us two small numbers 
a s (M) and v 2 that we can exploit to understand heavy quarkonium physics. 
In order to exploit the smallness of a s (M) , we need to solve the problem of sep- 
arating the scale M from the smaller momentum scales Mv, Mv 2 , and Aqcd 
that involve nonperturbative physics. In order to exploit the smallness of v 2 , 
we need some way of organizing nonperturbative effects according to how theii 
scale with v. These problems were solved by Peter Lepage and collaborator!! 
by introducing an effective field theory called nonrelativistic QCD (NRQCD). 

NRQCD is a formulation of QCD in which heavy quarks are treated as non- 
relativistic particles. Instead of describing the heavy quark and antiquark by a 
single 4-component Dirac field 'J, they are described by separate 2-component 
Pauli fields i/> and x ■ The lagrangian for NRQCD has the form 

•CnRQCD = Aight + Co + C% + • • • , (6) 

where Aight is the usual QCD lagrangian for the gluons and light quarks and 
antiquarks, Cq is the minimal NRQCD lagrangian, 

^ = ^{iDo + ^)^ + x'(iDo-^)x, (7) 

and Li, A, ■ ■ ■ are a series of "improvement terms" that can be added to 
make NRQCD reproduce the physics of full QCD to higher and higher order 
in v 2 . The minimal QCD lagrangian Cq is invariant under heavy-quark spin 
symmetry, which mixes heavy quarks with spin up and spin down. The v 2 - 
improvement term £2 includes four terms: a (D 2 ) 2 term that includes the 
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Table 2: The QCD coupling constant at the momentum scales for quarkonium systems 





cc 


bb 


tt 


a a (M) 


0.35 


0.22 


0.11 


ct s (Mv) 


0.52 


0.35 


0.16 


a a (Mv 2 ) 


£ 1 


£ 1 


0.35 



relativistic correction to the kinetic energy, a D • E term, and D x E • er and 
B • er terms that break the spin symmetry. The coefficients of these four terms 
depend only on the momentum scale M and they can therefore be calculated 
as perturbation series in ot s (M). 

The NRQCD collaboration has calculated the spectrum for bottomonium 
and charmonium using Monte Carlo simulations of lattice NRQCDel These cal- 
culations demonstrate convincingly that NRQCD provides an effective frame- 
work for describing heavy quarkonium physics. Using the lagrangian £ii g ht+A)> 
the splittings between spin-averaged energy levels are reproduced to an ac- 
curacy of about 30% for charmonium and 10% for bottomonium, which are 
consistent with errors of relative order v 2 . Because this lagrangian has spin 
symmetry, it gives no spin splittings. After adding the v 2 -improvement term 
£2, the accuracy of the spin-averaged splittings is improved to 10% for charmo- 
nium and 1% for bottomonium, consistent with errors of order u 4 . This term 
also gives the correct spin splittings up to errors of order v 2 . Higher accuracy 
could presumably be achieved by adding the v -improvement term £4. 

NRQCD provides a simple solution to the problem of separating the scales 
M from the lower momentum scales Mv, Mv 2 , and Aqcd- The improvement 
terms in the NRQCD lagrangian have coefficients that must be tuned to re- 
produce the physics of heavy quarks in full QCD. These terms compensate 
for the fact that NRQCD does not treat correctly the physics of relativistic 
heavy quarks whose momenta are of order M. Their coefficients therefore de- 
pend only on the momentum scale M. Since they depend only on the physics 
at length scales of order 1/M, we call them "short-distance coefficients." In 
NRQCD, all effects of the scale M are taken into account through the coeffi- 
cients of the improvement terms in the NRQCD lagrangian. Thus the NRQCD 
lagrangian itself provides us with the desired separation of scales. The effects 
of the scale M are all contained in the short-distance coefficients. 

NRQCD also provides a solution to the problem of organizing nonpertur- 
bative effects in such a way as to exploit the smallness of v 2 U It leads to a simple 
picture of the structure of a charmonium state, which is described most easily 
in terms of the Fock state expansion in Coulomb gauge. This gauge makes 
the dynamics of a nonrelativistic bound state particularly transparent, and it 
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allows for a sensible Fock state expansion since there are no negative norm 
states. For example, the Fock state decomposition for the J/tp in Coulomb 
gauge has the form 



\JH>) = E^Jcc^)) + E<i1^) 9 l cg 8( 3 ^)+5>+--- (8) 



The probability of each Fock state scales in a definite way with v. Only the 
cci( 3 <Si) Fock state has a probability of order 1. All the higher Fock states 
have probabilities suppressed by powers of v. The higher Fock state with the 
greatest probability is ccg( 3 Pj) + g, whose probability is of order v 2 . The 
next most important Fock state is ccg( So) + 9, whose probability is of order 
v 3 . One of the next most important Fock states is ccs( 3 Si) + g, whose prob- 
ability is of order v . By exploiting this Fock space structure, one can derive 
velocity-scaling rules that determine how matrix elements of local operators in 
a quarkonium state scale with v. If nonperturbative effects can be organized 
into matrix elements, their relative magnitudes can be estimated using the 
velocity-scaling rules. 

3.3 NRQCD Factorization Formula for Decays 

The NRQCD factorization approach to annihilation decays involves separating 
the scale M from the lower momentum scales Mv, Mv 2 , and Aqcdj and then 
exploiting the fact that quarkonium physics involves two small numbers, a s (M) 
and v 2 . The smallness of a s {M) is exploited by calculating the effects of the 
scale M as perturbation series in a s (M). The smallness of v 2 is exploited by 
using the velocity-scaling rules to estimate the magnitudes of nonperturbative 
matrix elements involving scales of order Mv and smaller. 

We proceed to outline the derivations of the factorization formula for an- 
nihilation decay rates. The inclusive rate for the decay of a quarkonium state 
H into light hadrons can be expressed in the following schematic form: 



where iPqq is the wavefunction for H to consist of a Q and Q and TQQ^ hadrons 
is the T-matrix element for producing a particular final state consisting of light 
hadrons. The sum over all final-state hadrons in ^ can be replaced by a sum 
over all final-state partons. This simply amounts to a change of basis for color- 
singlet final states. The final-state partons can be separated into hard partons 



cc 



ccg 




QQ hadrons 



(9) 
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with momenta of order M and soft partons with much smaller momenta. The 
resulting expression for the decay rate has the form 



T{H) = 2M^ E E EKo® T W-h ar d +sof t • (10) 

QQ hard soft 

The soft gluons in the final state can be emitted from the Q or Q, from the 
hard partons, or from the soft partons in the final state. The T-matrix element 
also involves virtual soft gluons that connect the initial Q and Q, the final hard 
partons, and the final soft partons. Standard factorization methods of pertur- 
bative QCD can be used to show that for every pair of hard partons whose 
4-momenta are not collinear, there is a cancellation between the corrections 
from virtual soft gluons exchanged between the two partons and the interfer- 
ence terms between real soft gluons emitted by the two partons. A similar 
cancellation occurs between the corrections from virtual soft gluons exchanged 
between a hard parton and the QQ pair and the interference terms from real 
soft gluons emitted by the hard parton and the QQ pair. The only soft gluons 
that survive after these cancellations are the virtual gluons exchanged between 
the Q and Q and the real gluons emitted from the Q and Q. But the effects of 
these gluons can be absorbed into the wavefunction ipQQ +so(t f° r H to consist 
of a QQ pair plus soft gluons. The resulting expression for the decay rate has 
the form 

r w = iEEE Ko+sott ® T QQ ^ haid \ 2 . (ii) 

QQ hard soft 

The sum over QQ states includes a sum over their color (1 or 8) and angular- 
momentum ( 2S+1 Lj) quantum numbers, which we denote collectively by n, 
and an integral over the magnitude of their relative momentum q. Because of 
the wavefunction factor, the integral over q has support only for q of order Mv. 
All dependence on the scale Mv can be removed from the T-matrix element 
by expanding T in powers of q. Denoting the expansion coefficients by T and 
ssociating the factors of q and the integral over q with the wavefunction factor, 
we obtain 

r (#) = E ( E I^Hhardl) ( fq L E KoM+softf) • (I 2 ) 

H n \hard / V 1 soft / 

The sum over n includes all color and angular-momentum states of a QQ pair. 
Each term in the sum is the product of a short-distance factor that involves only 
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the scale M and a long-distance factor that involves only scales of order Mv 
or smaller. The short-distance factor is proportional to the annihilation rate 
at threshold for a QQ pair in the state n. The long-distance factor measures 
the probability of finding the QQ pair at the same point in the state n in the 
quarkonium state H. 

The factorization formula ( |l2| ) is simply a statement about the separation 
of scales in QCD, so it does not necessarily require NRQCD. However NRQCD 
provides a convenient prescription for carrying out this separation of scales. 
The long-distance factor in ( |l2[) can be expressed as the expectation value in 
the quarkonium state H of a local gauge-invariant operator O n in NRQCD. The 
short-distance factor is proportional to the imaginary part of the coefficient of 
the operator O n in the NRQCD lagrangian. Thus ( |l2| ) can be written 

T W = ^ E ^QQ[n]{H\O n \H). (13) 
n 

This is the NRQCD factorization formula for the decay rate. Since the coeffi- 
cient r involves only the scale M, it can be calculated as a perturbation series 
in a s (M). The matrix elements (H\O n \H) can in principle be calculated non- 
perturbatively using Monte Carlo simulations of lattice NRQCDLD Thus the 
NRQCD factorization formula can be used to calculate the decay rate from 
first principles. 

The NRQCD factorization formula ([l3|) contains infinitely many terms. 
The utility of the formula relies on the fact that NRQCD also provides a 
way of exploiting the smallness of v 2 . Each of the matrix elements scales 
with a definite power of v. In the case of J/tp, the largest matrix element is 
(J/ip\Oi( 3 S±)\ J/tp), which is proportional to the square of the wavefunction 
at the origin and scales like v 3 . The next most important matrix elements 
are the expectation values of OsCSq), which scales like v 6 , and of 0$( 3 Si) 
and Os( 3 Pj), which scale like v 7 . The relative magnitudes of the terms in 
the factorization formula ( |l3|) are determined by the order in v of the matrix 
elements and by the order in a s of the short-distance coefficients. 

3-4 NRQCD Factorization Formula for Production 

The NRQCD factorization approach to quarkonium production is based on the 
fact that a sufficiently inclusive cross section satisfies a factorization formula 
analogous to ( |l^ ) . All effects involving momentum scales of order M and larger 
are contained in short-distance factors that involve the T-matrix element for 
producing a QQ pair in the state n with given 4-momentum P. All effects 
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involving momentum scales of order Mv and smaller are contained in long- 
distance factors that involve the amplitude for a QQ pair in the state n to form 
the quarkonium state H plus soft partons. The effective field theory NRQCD 
provides a convenient prescription for carrying out this separation of scales. It 
can also be used to exploit the small number v 2 by using velocity-scaling rules 
to determine how the long-distance factors scale with v. 

The factorization formula for the differential cross section for producing a 
quarkonium state H with 4-momentum P has the form 

da(H(P)) = Y, dd QQln,P](°n), (14) 

n 

where the sum over n includes all color and angular-momentum states of the 
QQ pair. The short-distance factor a is proportional to the cross section for 
producing a QQ pair at threshold in the state n with total 4-momentum P. 
Since they involve only momentum scales of order M and larger, the short- 
distance factors can be calculated as perturbation series in a s (M). The long- 
distance factor (O^), which can be expressed as an NRQCD matrix element, is 
proportional to the probability that a QQ pair produced at a point in the state 
n will form the quarkonium state H plus soft hadrons, whose energies in the 
H rest frame are of order Mv or smaller. The derivation of the factorization 
formula breaks down if the 4-momentum of the quarkonium is collinear with 
that of any hadrons in the initial state. It therefore does not apply in the 
diffractive region. 

The NRQCD factorization formula ([l3|) contains infinitely many terms. 
The relative importance of these terms depends on several factors. The rela- 
tive magnitudes of the matrix elements depends on how they scale with v. The 
relative magnitudes of the short-distance coefficients depends on their order in 
a s , and also on how they scale with dimensionless ratios of kinematic variables, 
such as m c /pT in the case of production at large transverse momentum. In 
practice, there are several matrix elements of phenomenological importance for 
any given quarkonium state. Unfortunately we have no effective nonpertur- 
bative methods for calculating these matrix elements, with a few exceptions. 
The exceptions are matrix elements like (Of^ ( 3 Si)) that can be related to the 
wavefunction of the dominant cc\ Fock state. The remaining matrix elements 
must be treated as phenomenological parameters. The predictive power of the 
factorization formula resides in the fact that the same matrix elements must 
describe inclusive quarkonium production in all high energy processes. 

The factorization formula ( filf ) is simply a statement about the separation 
of scales in QCD. One can ignore the prediction of NRQCD for the relative 
magnitudes of the matrix elements, and simply take the formula as a model- 
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independent framework for analyzing quarkonium production. Any model that 
is consistent with QCD at short distances must be expressible in this form for 
some choice of the NRQCD matrix elements. In particular, the color-singlet 
and color evaporation models can be expressed as specific assumptions about 
the NRQCD matrix elements. The color-singlet model reduces to the assump- 
tion that the only matrix element that is nonzero is the color-singlet matrix 
element whose angular-momentum quantum numbers correspond to those of 
the cc pair in the dominant Fock state. This matrix element is (O^" ( 3 Si)) 
for the J/4> and (p\ cj { 3 Pj)) for the Xcj- The color evaporation model corre- 
sponds to the assumption that the S-wave matrix elements dominate and that 
they are equal up to simple color and spin factors: 

(0f( 3 Si)) = (OffSo)) = 6<0f( 3 Si)) = 6(Of(^ )). (15) 

The P-wave matrix elements are all suppressed by a factor of (m 2 D — m 2 )/m 2 , 
which is of order v 2 . NRQCD predicts a much more intricate pattern for the 
relative magnitudes of the matrix elements and this pattern depends on the 
J quantum numbers of the quarkonium state H . The matrix elements that 
scale with the leading power of v include both the color-singlet model matrix 
element and one of the color evaporation model matrix elements in ( |l5| ) . For 
the J/ip, the only leading matrix element is (Of^ ( 3 Si)), which scales like v 3 . 
For the Xcj, the leading matrix elements are {0* cj ( 3 Pj)) and (Og CJ ( 3 Si)), 
both of which scale like v 5 . 

As a phenomenological framework for calculating quarkonium production, 
the NRQCD factorization approach is much less economical than the color 
evaporation model. For each quarkonium state, there are several matrix el- 
ements of phenomenological importance. For J/ijj production, the largest 
matrix element is (0(^ \ 3 Si)} , which scales like v 3 . Its value can be de- 
termined from decays of the J/ip. The next most important matrix elements 
are (O^ i 1 S )) , (0^{ 3 Si)), and (0^( s Pj)), which scale like v 6 , v\ and 
v 7 , respectively. These three color-octet matrix elements must be treated as 
independent phenomenological parameters, in contrast with the single phe- 
nomenological parameter /m in the color evaporation model. 

The NRQCD factorization approach differs dramatically from the color 
evaporation model in its predictions for the dependence of the cross section on 
the polarization of the quarkonium state. The color evaporation model pre- 
dicts no dependence on the polarization. The NRQCD factorization approach 
predicts a nontrivial dependence, because the NRQCD matrix elements (O^) 
depend on the polarization of the quarkonium state H. 
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4 Applications to Charmonium Production 

The NRQCD factorization approach has been applied to charmonium produc- 
tion in almost all possible high energy processes. Some of these applications 
are described in two recent reviews™ I will focus on one specific application for 
which the predictions of the NRQCD factorization approach differ significantly 
from both the color-singlet and the color evaporation models. 

4-1 J/ip at Large Transverse Momentum 

The NRQCD factorization formula leads to a dramatic prediction for the pro- 
duction of prompt charmonium at large transverse momentum in pp collisions. 
The prediction is that at sufficiently large pr, most of the J/ip's will be trans- 
versely polarized. This prediction follows from several simple steps, which I 
proceed to discuss in detail. 

The first step in the argument involves the separation of the scale px from 
the lower momentum scales of order m c ancLsmaller. According to standard 
factorization theorems of perturbative QCD[L3 the inclusive production cross 
section for any hadron at large pr is dominated by fragmentation, which means 
that it is produced by the hadronization of a single high-py parton. As pointed 
out by Braaten and Yuanlid this factorization theorem applies to charmonium 
states as well as to light hadrons provided that pt m c V3 In the case of the 
J/tp, the factorization formula is 

da(pp J/V(P) + X) = V f dzda(pp -> i(P/z) + X) D^ J/4 ,(z) , (16) 

i J o 

where da(pp — > i + X) is the inclusive cross section for producing a parton 
of type i with larger transverse momentum pt/z. This cross section can be 
calculated in terms of the parton distributions for p and p and parton cross 
sections da that involve only the scale pr/z and can therefore be calculated 
as power series in q: s (pt/z). The fragmentation function D^j/^ gives the 
probability that the hadronization of the parton i produces a J/ip carrying a 
fraction z of the parton momentum. All effects of momentum scales smaller 
than px, including the effects of the charm quark mass and nonperturbative 
effects involved in the formation of the bound state J/ip, are contained within 
the fragmentation functions. The largest fragmentation probabilities are those 
for c, c, and the gluon. The sum over partons i in ([R3]) is completely dominate 
by the gluon, because the cross section for producing gluons is so much larger 
than that for producing charm quarks. The corrections to the factorization for- 
mula (|l6|) from nonfragmentation processes are suppressed asymptotically by 
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powers of rri^/p^. Thus the factorization of the scale pr leads to the conclusion 
that J/ip production at large px is dominated by gluon fragmentation. 

The second step in the argument involves the separation of the scale m c 
from the lower momentum scales of order m c v and smaller. The NRQCD 
factorization formula ([l4]) can be used to express the gluon fragmentation 
function in the form 

D„^ m {z) = J2^cc[n)(z) (Of*) , (17) 



where d g ^ ce ^ is the fragmentation function for producing a cc in the state 
n carrying a fraction z of the momentum of the gluon. This fragmentation 
function involves only momenta of order m c , and it can therefore be calculated 
as a power series in a s (m c ). All nonperturbative effects involving the binding 
of the cc pair to form a J/tp are factored into the NRQCD matrix elements. 

According to the color-singlet model, the gluon fragmentation function 
( |l7| ) should be dominated by the cci( 3 Si) term. The leading contribution 
to the short-distance coefficient for this term is of order a 3 and comes from 
the parton process g* — > ccgg. Keeping only this term in the fragmentation 
function (|l7|), the cross section predicted by ( |l6| ) is about a factor of 30 below 
recent CDF data on prompt J /if) production at the Tevatron. It is this result 
that spelled the final demise of the color-singlet model. 

According to the color evaporation model, the gluon fragmentation func- 
tion (|l7]) should be dominated by the S-wave matrix element with the largest 
short-distance coefficient. We will see that the NRQCD factorization approach 
leads to the same conclusion. As first pointed out by Braaten and Fleming|i3, 
the term that dominates is the ccs( 3 Si) term, which corresponds to the for- 
mation of a J/tp from a cc pair that is created in a color-octet 3 Si state. At 
leading order in a s , this term in the fragmentation function is 



V^W « ^S^ 1 - ^^(^i))- (18) 



The relative importance of the various terms in the fragmentation function (17) 
depends on several different factors. We will compare each of these factors for 
the cci( 3 Si) term and the ccg( 3 5'i) term. The magnitude of the NRQCD ma- 
trix element is determined by its order in v, which is v for (0(^( 3 Si)) and 
v 7 for (Og ( 3 Si)). The magnitude of the short-distance coefficient d is de- 
termined by its order in a s . The coefficient of the cci( 3 <Si) term comes from 
the parton process g* — > ccgg and is of order a 3 , while the cci( 3 S'i) term has 
a much larger coefficient of order a s coming from the parton process g* — > cc. 
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However the relative importance of the terms also depends on the shape of the 
z-distribution of d. The reason is that in the expression ( [l6| ) for the cross sec- 
tion, the fragmentation function is folded with the cross section for producing 
a gluon, which is a rapidly falling function of the transverse momentum pr/z 
of the gluon. Since da jdp\ falls like the inverse fourth power of the gluon's 
transverse momentum, the cross section is effectively weighted by z 4 . This dra- 
matically suppresses the cci( 3 Si) term, which has a rather soft z distribution 
with (z) « 0.4. In contrast, there is very little suppression of the ccs( 3 Si) term 
in Jig), since it is sharply peaked near z = 1. Putting all the factors together, 
the cci( 3 5'i) term scales like (z A )a 3 v 3 while the ccs( 3 Si) term scales like a s v 7 . 
The suppression of the color-singlet term by {z A )a 2 s overwhelms the suppres- 
sion of the color-octet term by v 4 . Thus the factorization of the scale m c leads 
to the conclusion that the gluon fragmentation function is dominated by the 
ccg( 3 Si) term. This term is singled out because it's short-distance coefficient 
is the lowest order in a s and the most sharply peaked near z = 1. 

Unfortunately, we cannot use the formula ( Jl8| ) to predict the cross section 
at large pr, because the matrix element (0$^( 3 Si)) has not yet been deter- 
mined accurately from other production processes. However, the shape of the 
cross section da/dpT obtained by inserting ( |l8| ) into ( |l7| ) is in agreement with 
the Tevatron data at large pt- The value of the matrix element can therefore 
be determined by fitting the data. The resulting value of (O^^ ( 3 Si)) is con- 
sistent with suppression by a factor of v relative to the color-singlet model 
matrix element (0( /4 '( 3 Si)) . 

The color evaporation model gives a prediction for the J /if) cross section 
at large pr that is similar to that of the NRQCD factorization approach, be- 
cause the color evaporation model matrix elements in ( p"5| ) include (0g ( 3 Si)}. 
However the NRQCD factorization approach differs dramatically in its predic- 
tions for the dependence of the cross section on the polarization of the J /if). 
The color evaporation model predicts the absence of any spin asymmetry. As 
pointed out by Cho and WiseE 3 ], the NRQCD factorization approach predicts 
that J/i/j's should be predominantly transversely polarized at sufficiently large 
Pt- The leading term (|lgj ) in the gluon fragmentation function contributes 
only to the production of J/^'s that are transversely polarized. The radia- 
tive corrections to the fragmentation function were examined by Beneke and 
RothsteinLj, and they concluded that the spin alignment at large pr remains 
greater than 90%. The largest corrections to the spin alignment at values of 
Pt that have been measured at the Tevatron come from nonfragmentation con- 
tributions that are suppressed by m 2 c /p T relative to the terms in Jl6|). These 

Sntributions have been calculated by Beneke and Kraemer and by Leibovich 
. Their prediction is that the spin allignment should be almost zero at small 
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Pt, that it should begin to turn on at a px of 3 to 5 GeV, and that it should 
reach 60% to 80% of its maximum possible value for pt of 20 GeV. An experi- 
mental measurement of the spin alignment in agreement with these predictions 
would constitute a dramatic triumph for the NRQCD factorization approach. 
It would also spell the final demise of the color evaporation model. 



5 Conclusions 

The NRQCD factorization approach has led to dramatic progress in our under- 
standing of inclusive quarkonium production. Within this approach, the roles 
played by color and spin in the production process are somewhere intermedi- 
ate between the extreme assumptions of the color-singlet and color evaporation 
models. In some cases, the NRQCD factorization approach gives predictions 
that differ dramatically from those of both the color-singlet model and the 
color evaporation model. A particularly important example is the production 
of J/ip at large pt in pp colliders, where a measurement of the spin asym- 
metry will provide a stringent test of our present understanding of inclusive 
quarkonium production. 
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